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Cd ' Abstract 

Given a non-rational real space curve and a tolerance e > 0, we present an al- 
gorithm to approximately parametrize the curve. The algorithm checks whether 
a planar projection of the space curve is e-rational and, in the affirmative case, 
generates a planar parametrization that is lifted to an space parametrization. 
This output rational space curve is of the same degree as the input curve, both 
have the same structure at infinity, and the Hausdorff distance between them is 
always finite. 
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Introduction 

The development of approximate algorithms for algebraic and geometric problems is 
an active research area (see e.g. [21] and [23]) that focuses on different problems as, for 
instance, the computation of gcds (see [1], [10], [16]), the factorization of polynomials 
([6], [12], [IS]), the implicicitization of surfaces ([7], [9]), the parametrization of curves 
and surfaces ([2], [11], [H], [17], [18], [20]), etc. These approximate algorithms are 
applicable by themselves since they face symbolic computation to real world problems. 
Moreover, those of geometric nature are of special interest in the field of CAGD. For 
instance, providing parametric representations of algebraic geometric objects helps in 
some CAGD constructions as surface-surface intersection or computation of planar 
sections (see e.g. Example 5.3. in [20]). 



In this context, when using the term "approximate", there is certain risk of am- 
biguity since it can have a double meaning (see e.g. [21]). Let us clarify what we 
mean in our case. Here, as in our previous papers [T7], [IB], [20], with the term "ap- 
proximate" we do not mean "numerical" but something different (see introduction to 
[20] for further details): our input is the perturbation of an unknown input; once the 
perturbation is received we treat it exactly to provide an output that is close (in this 
sense it is "approximate"), under certain distance, to the theoretical output of the 
unperturbed unknown input. 

More precisely, the problem in this paper is as follows. Let C* be a rational real 
space curve defined as the complex-zero set of a finite set A^ C M[x,?/,2] of real 
polynomials; in practice card(AI) = 2. Nevertheless, instead of getting A4 as input 
of our problem, we get a new finite subset J-" C M[x, y, z] (which is a perturbation of 
Ai) of real polynomials that defines a new curve C, obviously different to C*. Since 
the genus of a curve is unstable under perturbations, the input curve C will have 
positive genus and hence it will not be parametrizable by rational functions. Ideally, 
the problem would consists in finding the initial curve C* or, even better, a rational 
parametrization of it. However, this goal is unrealistic. Instead, one might require to 
find a rational parametrization of the closest rational curve to C under certain distance; 
say, under the Hausdorff distance. Nevertheless, in this paper we deal with a weaker 
statement of the problem. Namely, finding a rational parametrization of one rational 
space curve being close (in comparison to a given tolerance e > 0) to C under the 
Hausdorff distance. Our statement, although may not yield to the best output rational 
curve, generates one good answer. This can be seen as a first step for the harder, more 
general, and theoretical problem of finding the best (in the sense of closest) rational 
curve being, our solution, meanwhile ready to be used in applications. 

In [IT], the authors show how to solve the problem for the particular case of 
e-monomial plane curves (i.e. plane curves having an e-singularity of maximum e- 
multiplicity); see also [IB] for the case of surfaces. Later, in [20], the problem was 
solved for the more general case of e-rational plane curves. The current paper is, there- 



fore, the natural continuation of this research since it deals with the next step, namely 
the case of space curves. 

In the unperturbed case, the problem can be solved by birationally projecting the 
space curve on a plane, checking the genus of the projected curve and, in case of genus 
zero, parametrizing the plane curve to afterwards inverting the parametrization to a 
rational parametrization of the input curve (see e.g. [22] for further details). Now, 
the situation is more complicated. More precisely the strategy (see box below) is as 
follows. 
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We get C, 
and close to C. 
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Applying Chinese Re- 
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We get a plane curve 
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We compute the 
e-genus of itz{C); 
say it is zero. 



Algorithm in [20] out- 
puts a rational curve V, 
parametrized as Q{t) = 
Cii/i§), and close to 

TTziC). 



We assume some conditions on the original space curve C (see Section [T|) such that 
when it is projected onto a plane we get a curve satisfying the hypotheses required by 
the algorithm in [20]; let us denote by 7r^(C) the projected curve, so we are assuming 
w.l.o.g. in this explanation that projection has been performed on the plane z = 0. 
Then, the algorithm in [20] determines whether iTziC) is e-rational, where e is a fixed 
given tolerance. If iTziC) is not e-rational one may try a different projection, but 
here we simply ask the algorithm to terminate since, although in some examples this 
seems to work, we do not have any theoretical argumentation to ensure when that 
projection exists. Otherwise, the algorithm in [20] goes ahead and computes a rational 
parametrization Q{t) of a plane rational curve T>. The last step consists in lifting V 
to a rational space curve C being close to C. For this purpose, we first realize that a 
sufficient condition for the finite Hausdorff distance requirement, between both curves, 
is given by the structure at infinity of the input curve C. Taking into account this fact, 
and using a Chinese-remainder type interpolation, we get a rational parametrization 



of C. As a consequence of this process, we get a rational space curve C of the same 
degree as C, having the same structure at infinity as C, and such that the Hausdorff 
distance between C and C is finite. 

The structure of the paper is as follows. In Section [H we introduce the notation 
that will be used throughout the paper as well as the general assumptions. Moreover, 
we comment on the reasons for the inclusion of these assumptions, we discuss how 
to check them algorithmically, and we show that they (the assumptions) are general 
enough. Section [2] is devoted to the projected curve iTz{C) and, more precisely, to prove 
that under the general assumptions imposed in[T]7r2(C) satisfies all requirements in the 
algorithm in [20]. Section [3] focuses on how to lift the rational plane curve V (generated 
by applying the algorithm in [20] to nz{C)) to the curve C such that both curves, C 
and C, have the same structure at infinity; note that vr^ is a birational map between 
C and vr^(C), but we are lifting V ^ 7r^(C). In Section |4] we summarize these ideas to 
derive an algorithm that is illustrated by two examples. In Section [5l we prove that the 
Hausdorff distance between the input and output curves, of our algorithm, is always 
finite. For this purpose, we briefly study the asymptotes of space curves. Finally, in 
Section El we approach the study of the Hausdorff distance empirically analyzing the 
examples in Section HJ 

1 General Assumptions and Notation 

We consider a computable subfield K of M, as well as its algebraic closure F; in practice, 
we may think that IK = Q. We denote by F^ and F'^ the affine plane an affine space 
over F, respectively. Similarly, we denote by P^(F) and P'^(F) the projective plane and 
projective space over F, respectively. Furthermore, if ^ C F^ (similarly if ^ C F^) 
we denote by A* C F^ its Zariski closure, and by A'^ C P'^(F) the projective closure 
of A*. We will consider {x,y,z) as affine coordinates and {x : y : z : w) as projective 
coordinates. Also, for A as above, we denote by A°° the intersection of A^ with the 
projective plane of equation w = 0. In addition, for every polynomial H G ]K[x, y, z] 
we denote by H^{x, y, z, w) the homogenization of H. 

Our method will be based on the projection of the space curve on a plane. Without 
loss of generality (see below) we will consider that 2; = is the projection plane. So we 
introduce the map 

TT, -.W^ ^¥^,{x,y,z) ^{x,y), 

as well as 

^h . p3(^p^) \ {(0 : : 1 : 0)} ^ P^(F), {x : y : z : w) ^ {x : y : w). 

Our main object of study will be an irreducible (over F) affine real (non-planar) space 
curve C C F^. Although, in practice, in most cases, C will be given by two generators, 
we present the results for the general case where a finite set of generators is provided. 



Therefore, we assume that C is given as the zero-set (over F) of a finite set of real 
polynomials J-" = {Fi,...,Fs} C ]K[x,?/, 2], s > 2. e will be the tolerance and we 
assume that < e < 1. In addition, we assume the following: 

General Assumptions 

1. The cardinality of C°° is deg(C). 

2. Tiz : C -^ T^ziCy is birational and deg(C) = deg(7r2(C)*). 

3. (1 : : A : 0), (0 : 1 : /i : 0), (0 : : 1 : 0) ^ C°° for any A, /i G F. 

4. If (1 : A : /i : 0), (1 : A : /i* : 0) G C°° then fx = /i*. 

5. The coefficient of Fi in z^'^'^^^^'^' is a non-zero constant; where tdeg denotes the 
total degree of Fi. 

We briefly comment on the reasons for the inclusion of the above assumptions, 
and we describe how to check them algorithmically. The condition on irreducibility 
is natural since rational curves are irreducible varieties. In any case, one can always 
consider the irreducible decomposition of the input to apply the results to each of the 
irreducible components. The assumption on the reality of the curve is included because 
of the nature of the problem, but the theory can be similarly developed for the case of 
complex non-real curves. The exclusion of planar curves is to simplify the exposition. 
Note that this is not a loss of generality, since one can always apply the algorithm in 

m- 

Concerning the general assumptions, condition (1) will play a fundamental role in 
the error analysis, and it will be used to ensure that the Hausdorff distance between out- 
put and input is always finite. The birational requirement in condition (2) is introduced 
to reduce the problem to a plane curve after projection, and the degree fact will be used, 
in combination with (3) and (4), to ensure that 7rz{C)* has as many different points at 
infinity as its degree; this condition is required by the algorithm in [20]. Conditions (3) 
and (4) are related to the projection vr^. On one hand, (1 : : A : 0), (0 : 1 : /i : 0) ^ C°°, 
for any A, /i G F, ensures that (1 : : 0), (0 : 1 : 0) ^ 7rz(C)°^ which is a requirement for 
the algorithm in [20] to be applied to 7r^(C)*. On the other, (0 : : 1 : 0) ^ C°° guar- 
antees that vr^ is well defined on C^. In addition, conditions (3)- (4) ensure that vr^ is 
injective on C°°. Condition (5) is also introduced to guarantee that iiziC)* satisfies the 
hypotheses in [20] (see Theorem 12. 4p . Note that this property can always be achieved 
by means of a suitable orthogonal affine change of coordinates, and hence preserving 
distances. 

Taking into account that, in practice, C is expected to come from the perturbation of 
a rational space real curve, in general, all above conditions will hold. Nevertheless, let 
us discuss how to decide algorithmically whether a given input satisfies them. Checking 
the irreducibility of C can be approached by checking whether the corresponding ideal 



is prime (see, for instance, Section 4.5 in [8] or [5]). In order to check the reahty, one 
can apply cyhndrical algebraic decomposition techniques to decide the existence of real 
regular points (see e.g. [3]). The non-planarity of C can be deduced from a Grobner 
basis of J-". Let us now deal with the general assumptions. One can compute C^ by 
homogenizing a Grobner basis of J-", w.r.t. a graded order (see e.g. page 382 in [S]). 
The degree of C can de determined by counting the number of intersections of C with 
a generic plane; in fact a randomly chosen plane might be enough. So, (1) can also 
be checked. Also, (3) and (4) are checkable. Condition (2) can be analyzed by direct 
application of elimination theory techniques. Condition (5) is trivially checkable. 

In the previous description, we have considered that ^r = is the projection plane. 
Indeed, conditions (3)- (5) depend on this fact. So, if any of these conditions fails 
we might either consider a suitable orthogonal affine change of coordinates or choose 
another projection plane. Also, if (2) fails we need to find a different projection plane. 
We recall that, for almost every plane, the corresponding projection is birational over 
C and that for almost every plane the number of intersection points of the plane with 
C is deg(C). Therefore, the combination of these two facts with Lemma ll.ll ensures 
that condition (2) must be achieved by taking the projection plane randomly. 

Lemma 1.1. Let 11 C F^ 6e a plane such that card(C fl 11) = deg(C), let u be a (non- 
zero) parallel vector to U and non-parallel to the vectors in {P — Q\P,Q^Cr\Il, P ^ 
Q}, and let 11" be any plane orthogonal to u. Then, deg(7rn«(C)*) = deg(C), where vrn^ 
is the projection map from F^ onto 11". 

Proof. Let d = deg(C) and C nU = {Pi,...,Pd}, and let L be the line U fl H". 
By construction, {vrn"(-Pj)}.j=i,...,rf C 7rn«(C)* fl L. Since u is not parallel to Pi — Pj, 
with i 7^ j, then card({7rn"(-Pi)}i=i,...,d) = d. Therefore, deg(C) < deg(7rn«(C)*). Now, 
let L' be a line in 11" such that L' fl (7rn«(C)* \ 7iu^{C)) = and such that card(L' fl 
7rn"(C)*) = deg(7rn"(C)*); note that almost all lines in 11" satisfy this property. Let 
L' n 7rn"(C)* = {Qi, . . . , Qd'} and let 11' be the plane containing L' and being parallel 
to u; note that u is normal to 11", and hence L' is not parallel to u. Because of the 
construction n^l{Qi) flC 7^ and it is contained in 11'. Therefore, U^=i%l((5i) flC fill' 
has cardinality at least d', and hence deg(7rn"(C)*) = d' < deg(C). n 

2 The Projected Curve 

In this section, we analyze the basic properties of the projected curve 'Kz{C)* . In 
particular, we show that it satisfies the hypotheses in [20]. We recall that, since tt^ is 
birational and C is irreducible, iTziC)* is irreducible. We start with a technical lemma 
on Grobner bases 

Lemma 2.1. Let L C F &e a field and Gi, . . . ,Gm G h[xi, . . . ,Xn] be such that 
deg^^{Gi) = tdeg(Gi) > 0; where tdeg denotes the total degree. Let {Hi,...,Hr} 
be a Grobner basis of (Gi, . . . , Gm) w.r.t. the graded lex order with Xi < ■ ■ ■ < Xn. 
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Then, there exists i G {l,...,r} such that deg^^(ifj) = tdeg(ifj). Moreover, if the 
variety defined by {Gi, . . . , Gm} over F is not empty then deg^^{Hi) = tdeg{Hi) > 0. 

Proof. Let i = deg2.^(G'i). Since i = tdeg(G'i), because of the ordering, the leading 
term of Gi is x^. Now, by Exercise 5, page 78, [8], there exists i G {1, . . . ,r} such 
that the leading term of Fj divides x^. Finally, because of the ordering, tdeg(ifj) = 
deg^^(ifj). Now, if the variety of {Gi, . . . , Gm)} is empty, we assume that the Grobner 
basis is normal (this does not affect to the previous reasoning). By Theorem 8.4.3 in 
[2F| . Hi is not constant. So tdeg(ifj) > 0. D 

The next lemma shows how generalized resultants can be used to compute the 
projection. 

Lemma 2.2. Let 

where A is a new variable, and let F^ be the homogenization of F^{x,y,w,A) as a 
polynomial in ]K[A] [x, y, z] . Let 

m m' 

R = Res,(Fi,FA) = ^a,(x,2/)A^ S = Res,(F^F^) = ^ A(x,y,w;)A\ 

It holds that 

1. Tiz{C)* is the affine plane curve defined by gcd(ao, . . . , Om)? O'^d m = m! . 

2. If J^ is a Grobner basis, w.r.t. the graded lex order with x < y < z, then tTz{C)^ 
is the projective plane curve defined by gcd(/3o, • • • , Pm)- 

Proof. (1) We first prove that 7iz{C)* is the variety defined by {ao, . . . ,am}- In- 
deed, let (a, 6) G iTz{C). Then, there exists c G F such that P = {a,b,c) G C. 
So, Fi(P) = 0,Fa(P,A) = 0. Thus, R{a,b,A) = 0, and hence ao{a,b) = ■■■ = 
am{a, b) = 0. Conversely, let oq, • . • , am vanish at (a, b). Then R{a, 6, A) = 0. Now, 
since deg^{Fi) = tdeg(Fi), there exists c in the algebraic closure of F(A) such that 
Fi{a,b,c) = 0,F/^{a,b,c,A) = 0. Since c is a root of Fi{a,b,z) G ¥[z] then c G F. 
Therefore, from F^{a,b,c,A) = 0, we get that ^2(0, 6, c) = ■■■ = Fs{a,b,c) = 0. So, 
(a, b,c) E C and (a, b) G 7r^(C). 

Let a = gcd(ao, . . . , Qm) and a^ be such that Ofj = «,«. Let V and W be the 
varieties defined by a and {ao, . . . ,am}, respectively. Then, nz{C)* = V U W. Since 
C is irreducible and tt^ birational, we have that 7iz{C) is an irreducible curve. So, V 
is 1-dimensional and W is either empty or 0- dimensional. In any case, because of the 
irreducibility, W C V. So 7r^(C)* = V. 

Finally, let us see that m = m'. We assume w.l.o.g. that all Pi are non-zero. Since 
deg^(Fi) = deg^(F/*), by Lemma 4.3.1. in [23], we get that R{x,y,A) = S{x,y,l,A), 
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up to multiplication by a non-zero constant. Moreover, S is homogeneous as a poly- 
nomial in F[A][x,|/, w]. Thus, /?« are homogeneous (of th e same degree). Therefore, 
/3j(x, y, 1) does not vanish. Thus, m = deg^(-R) = deg^(S'(x, y, 1, A)) = deg^(S') = m'. 

(2) We first prove that w does not divide S. Let S = wM{x,y,w,A). Then, for 
all (a, 6) G F^, since deg(F/') = deg^(F/'), there exists c G F (in principle, c is in 
the algebraic closure of F(A) but, reasoning as above, we get that c G F) such that 
Fl'(a,b,c,0) = F^(a,6,c,0,A) = 0. Therefore, since c G F, i^'^(a, 6, c, 0) = 0,i = 
1, . . . , s; let us call p{a, b) the corresponding c associated to (a, b). Then, the infinitely 
many points {(1 : n : p{l,n) : 0)}„gN are included in the intersection of C'^ with the 
plane w = 0, which is a contradiction; note that since J-" is a Grobner basis w.r.t. a 
graded order then {F^ \ i = 1, . . . ,s} generates C^. 

From R{x, y, A) = S{x, y, 1, A) we get that aj{x, y) = /3j{x, y, 1). Therefore, since 
13 j is homogeneous, a^w^^ = /3j, for some rij G N. Moreover, since w does not divide 
S, there exists io ^ {0, . . . , m} such that a^^ = /3ig and gcd(afp, w) = 1. 

Let a = gcd(Q!o, . . . , «„,) and 7 = gcd(aQ, . . . , a^). We see that a^ = 7. Let 
ttj = aUi. Then, a^ = a'^a^. So, a^ divides 7. Conversely, let a^ = 'jai. Then, 
7(x, y, 1) divides a^{x,y,l) = ai. Therefore, 7(x,?/, 1) divides a. In addition, since 
a^ divides 7, a divides 7(x,y, 1). Hence, up to multiplication by non-zero constants, 
a = 7(x,y, 1). Therefore, since by construction w does not divide 7, we get that 
a'^ = 7. 

Finally, it remains to prove that 7 = gcd(/3o, . . . , /3m)- We know that 
gcd(/3o, . . . , /3m) = gcd(aow"°, . . . , a'^w'^"'). Let a = gcd(aow"°, . . . , a'^w'^"'). Clearly 
7 divides a. Conversely, a divides af^ (see above). Since gcd(afp,w) = 1, then 
gcd(a,w) = 1. Therefore, a must divide all aj. Hence, a divides 7. 

Summarizing gcd(ao, . . . , ctm)^ = gcd(aQ, . . . , a^) = gcd(aQw"°, . . . , a^w"™) = 
gcd(/3o,...,/3m). D 

Remark 2.3. We observe that in the proof of Lemma 12.21 from all the hypotheses 
imposed in Section [H we have only used the following: general assumption (5) is used 
in both (1) and (2). The fact that C has dimension 1 and that is irreducible is used in 
(1), jointly with the fact that vr^ is finite. Finally, in (2), we use that C'^ intersects the 
plane w = in finitely many points. 

We finish the section by stating the main properties of the projected curve. 
Theorem 2.4. It holds that 

1. 7r,^(C°°) = 7r,(C)°°. 

2. card(7r,(C)-)=deg(7r,(C)*). 

3. (l:0:0),(0:l:0)^7r,(C)°°. 



Proof. (1) Because of Lemma |2.H we can assume w.l.o.g. that {Fi, . . . ,Fs} is a 
Grobner basis w.r.t the graded lex order with x < y < z, and that deg^(Fi) = tdeg(Fi). 
Also, let Fa, F^, S, R, ai, /3i be as in Lemma [2l2| and a = gcd(ao, . . . , dm)- 

Note that C°° is the zero set in P^(F) of {F/'(x, y,z,0),. . ., F^{x, y, z, 0)}. So, since 
C°° is zero-dimensional, then gcd(F/^(x, y, z,0), . . . , F^{x, y, z, 0)) = 1. In addition, by 
Lemma [2^ vr2(C)°° is the zero set in P^(F) of a'^{x, y, 0). 

Now, let (a : 6 : c : 0) G C°^. Then, F/^(a,6,c,0) = F^(a,6,c,0) = 0. Therefore, 
S{a,b,0, A) = 0. Thus, /3j(a,6, 0) = 0. By Lemma |2.2[ we know that a^{x,y,w) = 
gcd{(3o, . . . ,(3m). So, (3i = a'^/Sj. Let us assume that a'^{a,b,0) ^ (i.e. that (a : 
b : 0) ^ 7r^(C)°°). By general assumption (3), a,b cannot be both zero. We assume 
w.l.o.g. that a = 1. Then, /3j(l,6, 0) = for all i. We now consider the polynomials 
Hi{y, z, w) = F!^{1, y, z, w) as well as the affine variety V defined by them. Note that, 
since (1 : 6 : c : 0) G C^, (6, c, 0) G X^ 7^ 0. Moreover, since C is irreducible, V is an 
irreducible curve. Furthermore, tdeg(iJi) = deg^(ifi) = deg2(-Fi) > 0. Furthermore, 
since C is not planar, T> is not a line perpendicular to the plane 2; = 0, so tTz is finite 
over V. Also, since C is not planar, V'^ intersection x = has only finitely many 
points. Furthermore, because of Lemma [2?n we can assume w.l.o.g. that {Hi, . . . ,Hs} 
is a Grobner basis w.r.t the graded lex order with w < y < z, and that deg^(ifi) = 
tdeg(i7i). Thus, V satisfies the hypotheses of Lemma [22] (see Remark l2.3p . Let H^ 
as in Lemma [2^ let T{y, w, A) = ReSz{Hi,HA), and let S, ai, a^ , /3j be as in the proof 
of Lemma 1221 Reasoning as in the proof Lemma (221 we get that degA(T) = degA(5') 
and that, if T = ^"q/^*^* then pi{y,w) = f3i{l,y,w). Moreover, by Lemma [221 we 
get that 7iz{'D)* is defined by p = gcd(po5 • • • ,Pm)', note that 7rz{'D)* is irreducible, and 
hence p is an irreducible polynomial. 

From pi{y,w) = (3i{l,y,w) = a'^ {1 , y , w) f3i{l , y , w) we get that a^{l,y,w) divides 
p{y,w). Also, since C is not planar, tTz^C)* is not a line, and hence a^{l,y,w) is not 
constant. Thus, since p is irreducible, we get that, up to multiplication by non-zero 
constants, p{y,w) = a^{l,y,w). Finally, from (6, c, 0) G V, we get that p{b,0) = 
a^{l,b,0) = 0, which is a contradiction. This proves that 7r^(C°°) C 7iz{C)°°. 

On the other hand, because of general assumptions (3) and (4) one has that vr^ 
is injective over C°°, and hence we get that card(7r^(C°°)) = card(C°°). Then, from 
general assumptions (1) and (2), we get that 7r^(C°°) = '7iz{C)°°. 

(2) Because of general assumptions (3) and (4), card(C°°) = caid^n^ {C°°)) and, by 
general assumptions (1) and (2), card(C°°) = deg(C) = deg(7r2(C)*). Now the proof 
ends by applying statement (1) in this theorem. 

(3) It follows from general assumption (3) and statement (1) in this theorem. n 

3 The Lifted Curve 

In Theorem 12.41 we have seen that, under the assumptions introduced in Section [H 
TTziC)* satisfies the hypotheses required by the parametrization algorithm in [20] • In 
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this situation, we apply algorithm in [20] to tTz^C)*. If iTziC)* is not e-rational, then 
we can not use 'n'ziC)* to parametrize C approximately by this method. However, it 
might be that there exists another projection such that the projected curve is e-rational 
and hence the method applicable to this other projection. Nevertheless, we have not 
researched in this direction leaving this as a future research line. So, let us suppose 
that 'JTz{C)* is e-rational, and let 

Q.. _ (pM pjty 

be the parametrization output by the algorithm in [20]. Let V be the rational plane 
curve parametrized by Q{t). We want to lift T> from F^ to a rational curve C in F^. 
For this purpose, in order to guarantee that the Hausdorff distance between C and C 
is finite (see Corollary 15. 5p . we will associate to P a rational curve C in F^ such that 
TTziC) = V, deg(C) = deg(C) and C°° = C°°. 

We know that iiziC)* and V have the same degree and the same structure at infinity 
(see Theorem 4.5. in [20]). Thus, by Theorem 12. 4[ P°° = 7r^(C°°). In addition, it also 
holds that deg(pj) < deg(g) = deg(P) (see proof of Lemma 4.2 in [20]). Moreover, 
by construction, gcd{pi,q) = 1 (see Step 10 in the algorithm in [20]). Furthermore, 
(pi(t) : P2(t) : q(t)) reaches all points in V°° (see proof of Theorem 4.5. in [20]). 
Therefore, since card(P~) = deg('D) (see Theorem [23] and note that V°° = 7r^(C~)), 
q(t) is square-free. Thus, if {^i, . . . , ^d} are the roots of q(t), 

D-=|fl:^^:0 

U pite) 

because of general assumption (4), for every i there exists a unique Xi ^^ such that 

iv pite) 

Note that, if {Gi, . . . , Gm} is a Grobner basis of {Fi, . . . , Fs} w.r.t. the graded lex 
order with x < y < z, then Xi is the root of 

gcdfGtfi,^,.,oy...,G::,fi,^,.,o 

Let Psit) be the interpolating polynomial such that ^3(^1) = Pi{^i)Xi: ior i = 1, . . . ,d; 
recall that q(t) is square-free. We then define C as the rational curve 

^f,. ^ fMt)_ Mti P3{t) \ 

^'> [ q{t) ' q{t) ' q{t) J ' 

Note that gcd{pi,p2,P3,q) = 1, g is square-free, deg(pi), deg(p2) < deg(g) and 
deg(p3) < deg(g). 

Taking into account the previous reasonings, we have the following theorem. 
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Theorem 3.1. The lifted curve C, defined as above, satisfies that 

1. C is rational. 

2. C^ = C°°. 

3. deg(C) = deg(C). 

We finish this section explaining how to compute C (i.e. the polynomial Psi^t)) 
without having to explicitly compute the roots of q{t). The idea is to adapt the 
Chinese Remainder interpolation techniques. Let {Gi, . . . ,Gm} be as above, and let 
?(^) ~ Y[i=i1ji^) be an irreducible factorization of q{t) over K. Now, for each qj we 
consider the field L = ]K(/i), where /i is the algebraic element over IK defined by qj{t), 
as well as the polynomial ring L[2;]. Let 

D,{z) = gcd (g1 f 1, ^, ., oV . . . , g:;, f 1, ^, z, 

where the gcd is taken in the Euclidean domain h[z]. Because of the previous reasoning, 
we know that Dj{z) can be expressed as 

where u &N. Let Cj(/x) be the polynomial expression oi bj{fi)aj{fi)~^pi{fi) as an element 
in L. On the other hand, for i ^ j, gcd{qi,qj) = 1. So, there exist Uij,Uj^i G K[t] such 
that Uijqi + Uj^iqj = 1. We introduce the following polynomial 

A(f) - i (^lit)UL2Ui,l{t)q^it) + ■ ■ ■ + Ci{t)]^Zlu,,,it)q,{t) if £> 1 

Then, we have the following result. 

Lemma 3.2. Pslt) is the remainder of the division of A{t) by q{t). 

Proof. If £ = 1 the result is trivial. Let i > 1 and let R{t), Q{t) be the remainder and 
quotient of the division of A(t) by q(t), respectively. Clearly deg(_R) < deg(g). Now, 
let ^i be a root of q{t); say w.l.o.g. that ^i is a root of qi(t). Then, by construction, 
ci(6) = XiPii^i)- Therefore, 

i?te) = Afe) - gte)Qte) = ^(6) = cite>2,ite)g2te) ■ ■■Ui,^{^^)qei^^)■ 
However, for k ^ 1, Uk,i{^i)qk{^i) = 1 - Mi,fcte)gi(6) = 1- So, R{^i) = Ci(^i) = 

XiPl{^i)=P3{^i)- □ 
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4 Algorithm and Examples 

In this section we collect all the ideas developed in the previous sections to derive the 
approximate parametrization algorithm, and we illustrate it by several examples. For 
this purpose, we assume that we are given a tolerance < e < 1 as well as an space 
curve C satisfying all the hypotheses imposed in Section [TJ Then the algorithm is as 
follows 

Algorithm 

1. Compute the defining polynomial of iiziC)* (apply e.g. Lemma I2l2|) . 

2. Apply to iTz{C)* the parametrization algorithm in [20]. If the plane curve is not 
e-rational exit returning no parametrization else let (^7^-, ^^%t) be the output 
parametrization. 

3. Apply Lemma [3.21 to determine Psit). 
4 Return (^^ ^^ S^) 



Remark 4.1. Note that, because of Theorem 13. ![ it holds that the rational curve 
output by the algorithm has the same degree and structure at infinity as the input 
curve. Also, as already mentioned in Section |3l if in step 2 we do not get e-rationality, 
it does not imply that under another projection one could not get an e-rational curve. 
However we can not guarantee theoretically when such a projection exists. 

In the following examples, the polynomial / defining tTz^C)* and the parametriza- 
tions Q{t) of V, and V{t) of C, are expressed with 10-digits floating point coefficients, 
but the executions have been performed with exact arithmetic; the precise data can be 
found in http://www2.uah.es/rsendra/datos.html. 

Example 4.2. Let C be the space curve defined by the polynomials 

718945312497 698623125001 .7...,.., ^iq8«^^78«qq 

Fi = X H y — 671015625 z + 13865578693 zy 

100 100 ^ ^ 

- 12118499950 zx + 24392628607x2/ - 18401807886?/^ - 1311877532 z^ 

431020499999 1675347948801 

F2 = X + y - 1609143200 z + 4365980240 z y 

25 100 

- 401217042 zx- 24936051360 y^ - 683547137 z^ + 24392628607 x^ 

and let e = Y^. One can check that C satisfies all the hypotheses imposed in Section 
[H Moreover, deg(C) = 4. The projected curve nz{C)* is defined by the polynomial (see 

Fig. a2D 
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/(x, y) = 5.192147942 ■ lO^^ xy - 2.214420657 ■ lO^^ y - 5.059350678 ■ lO^^ x - 
2.636990684 • 10^^ x^ - 3.506554787 ■ 10^^ ^2y ^ 2.001041491 ■ 10^^ yi _ 
1.375243688 ■ 10^^ yS ^ 1. 181135404 • 10^^ ^3 ^ 3.822854018 ■ 10^^ ^y2 _ 
2.315025392 ■ W° y^ - 2.990857566 ■ lO^^^-^s _ 1.221346211 ■ 10^^3,2^2 ^ 
3.915698981 ■ 10^^ ^3y _ 1.812915331 ■ 10^^ ^,4^ 

Note that deg(7r^(C)*) = 4. Applying the approximate parametrization algorithm 
for plane curves in [20] we get that iiziC)* is e-rational. Furthermore the algorithm 
outputs the parametrization Q{t) = i^^, ^7^) 

-0.4173571408 + 1.171283433 i - 0.8477221239 1^ - 0.1445883061 1^ + 0.2133409452 1^ 
-0.9059858774 + 1.956830479 1 - 0.6103552658 f^ - 1.494650450 1^ + t^ ' 

0.1828752070 1 + 0.6268800173 1^ - 1.028340444 1^ + 0.3822448988 1^ - 0.1884116000\ 
-0.9059858774 + 1.956830479 1 - 0.6103552658 1^ - 1.494650450 1^ + t^ ) 

It only remains to compute the numerator of the third component of the rational 
parametrization of the lifted curve; namely P3{t). Applying Lemma 13.21 we get the 
approximate parametrization (see Fig. 14.21 P(t) = (^^77^-, ^77^-, ^jt)) of the space curve 
C (i.e. the parametrization of C), where 

Psit) _ -1.067157288^3- 0.2783759249 -0.7182447737t+ 1.955832944^2 
~qit) ~ -0.9059858774 + 1.956830479 1 - 0.6103552658 1^ - 1.494650450 t'^ + t^' 

Example 4.3. Let C be the space curve defined by the polynomials 

Fi =20052827033x1/ + 2850904342^1/ - 7155364672za; — x 

- 7869010116^ + 1 3412651801 ^ _ 43;L02722226y2 + 1610946062^^ 

390188402999 
F2 = - 18330943984^2/ + 33857630124zx x - 56921602320z+ 

12611223036001 , , , 
— y - 1666085147601/2 + 57179742076^^ ^ 20052827033^^ 

and let e = ^. One can check that C satisfies all the hypotheses imposed in Section 
[TJ Moreover, deg(C) = 4. The projected curve 7r^(C)* is not e-rational, but 7iy{C) is. 
So we work with the projection on the plane y = 0. The defining polynomial of 7Ty{C)* 
is 

f{x, z) = 6.959832072 ■ 10^^ zx - 4.075715387 ■ lO^^ x - 6.207866771 ■ 10=^^ z + 
1.769623619 ■ 10^^ x'^ + 9.541705261 ■ 10^*^ x^ + 1.269145848 ■ 10^^ z^ + 
8.077561390 ■ lO^^x^z - 1.355904241 ■ lO^^x^z - 2.573514563 ■ 10^^ z^ + 
2.289865008 ■ 10^^ z^ - 3.292700550 ■ 10^^ z^x + 4.798217962 ■ 10^^ xz^ + 
3.090311649 • W^x^z^ - 2.981944666 ■ 10^^ x^ 

13 







Figure 1: Up left: plot of C; Up right: plot of C; Down left: plot of vr^(C); Down right: 
plot V. 

Note that deg(7ry(C)*) = 4. Applying the approximate parametrization algorithm 
for plane curves in [20] we get that iiyiC)* is e-rational. Furthermore the algorithm 

outputs Q(t) = (a|),^) 

p^(t) = - 1.304559082 ■ 10^^^ t - 2.995071314 ■ 10^^^ t^ - 2.096039950 ■ 10^^^- 
1.114733279 ■ 10^^^ t^ - 3.005548232 ■ 10^^^ t^ 



P3{t) =6.238758852 ■ 10^^^ t - 5.937809784 ■ 10^^^ t^ + 3.175932541 ■ 10^^^ 
- 3.618736499 ■ 10^^^ t^ - 2.083861701 ■ lO^^'^t^ 
q(t) =3.555348439 1^ + t^ + 4.622830832 1^ + 2.625458073 1 + 0.5529230644. 



It only remains to compute the numerator of the second component of the rational 
parametrization of the lifted curve; namely P2{t)- Applying Lemma 13.21 we get the 
approximate parametrization V{t) = (^^^, ^^, ^^) of the space curve C (i.e. the 
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Figure 2: Joint plot of C and C 



parametrization of C) where 



P2{t) 



705 



= = 8.923772403 ■ 10 
5.846606727- 10^°^ t^ 



t + 1.249180934- 10^°^ t^ 



2.137471224-10'"°+ 



705 



5 Error Analysis 

In this section, we prove that the Hausdorff distance between the input and output 
curves, of our algorithm, is always finite. For this purpose, we first need to develop 
some results on asymptotes of space curves. Afterwards we will analyze the distance. 
To start, we briefly recall the notion of Hausdorff distance; for further details we refer 
to [1] . In a metric space {X, d) , for 7^ 5 C X and a G X we define 

d{a,B) = mib(^B{d{a,b)}. 

Moreover, for A,B G X \^ we define 

Hd(A, B) = max{sup„g^{d(a, B)}, supfeg5{d(6. A)}}. 

By convention Hd(0,0) = and, for 7^ A C X, Hd(v4,0) = cx). The function Hj 
is called the Hausdorff distance induced by d. In our case, since we will be working in 
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(C'^,d) or (R'^,d), being d the usual unitary or Euclidean distance, we simplify the 
notation writing H(A,i?). 

Let £ be an space curve in C^; similarly if we consider the curve in C". The intuitive 
idea of asymptote is clear, but here we formalize it and state some results. Although 
these results might be part of the background on the theory of asymptotes, we have 
not been able to find a reference in the literature suitable for our needs. 

We say that a line C in C'^ is an asymptote of £ if there exists a sequence {PnjneN 
of points in £ such that lim„||P„|| = oo and lim„d(P„,£) = 0; where d denotes the 
usual unitary distance in C^ and || || the associated norm. In the following, we show 
how the tangents at the simple points at infinity of £ are related with the asymptotes. 
More precisely, we have the following lemma. In the sequel, if V is a projective variety 
in P^(F), we denote by Va the open set V fl {tf 7^ 0}. In addition, for A G C, A denotes 
its conjugate. 

Lemma 5.1. Let P = (a : b : c : 0) & £°° be simple and let T{P,£^) be the tangent 
line to £^ at P. IfT{P,£^) is not included in the plane w = 0, then T{P,£ 
asymptote of £. Moreover, {a,b,c) is a direction vector of the asymptote. 



^^ IS an 



Proof. Let {ifi, . . . , ifm} C F[x,|/,2,w] be homogeneous polynomials defining the 
ideal of £^'^. Let 

7ri(x, y, z, w) = -K-{P)x + -^{P)y + -^{P)z + -^{P)w. 
ox oy oz aw 

Then, T(P, £^) is the projective variety defined by {ni{x, y, z,w), . . . , 7Tm{x, y, z, w)} 
(see e.g. pp. 181 in [13]). Note that, since T{P,£'^) is not included in the plane w = 0, 
T{P,£'^)a is an affine line. Now, we consider a local parametrization V{t) = {x(t) : 
y(t) : z(t) : w(t)) oi £^ centered at P. Since P is simple, and its tangent is not included 
in w = 0, the multiplicity of intersection of £'^ and the plane w = at P is 1. Thus, 
w(t) can be expressed as w(t) = tu(t), where u(t) has order 0. Therefore, u{t)~^ is a 
power series and Vit) can be expressed as 

V{t) = (x(t) : y{t) : z{t) : t) 

where x(t) = x(t)u~^(t), y(t) = y(t)u~^(t), z{t) = z(t)u~^(t) are power series. More- 
over, since T(P, £^) is the tangent, the order of HiiVit)) has to be, at least, 2. Therefore, 
TTi{V{t)) can be expressed as 7ii(V(t)) = t^^v{t), where ii > 1 and v(t) has order 0. 

Now, let {tn} be a sequence of complex numbers converging to 0, and such that 
tn 4 0. Then, for all n, P„ = (4^, ^, ^) e £. Moreover, lim„||P„|| = 00. Let H^ 

t^ri Cn <^n _^ 

be the affine plane defined by 7rj(x, y, z, 1). We prove that, for all 2, lim„d(P„, IIj) = 0. 
From where, one deduces that lim„d(P„, T(P, i''')^) = 0, and hence that T(P, £^'*)a is 
an asymptote. 
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Indeed, 



d(P„,n, 



ki(P„,i)| 



\ttMtn) 



dHj 
dx 



in^ip),^ip) 



dH,. 



dy 



dz 



dx 



dH, 



in^ip),^ip) 



dy 



dz 



Since the denominator is not zero, because T(P, S^) is not included inw = 0, and since 
ii — 1 > and order of v is 0, we get that hm„d(P„, IIj) = 0. 

Finally, we see that u = (a, 6, c) is a direction vector of the asymptote. First we 
observe that u is not zero, since P G P^(F). Now, by Euler equality we have that 



dx 
Substituting by P, we get 



-X + -Ti — y 



dy 



dHi dHi 

-z + — — w 



dz 



dw 



deg{Hi)Hi. 



dHj , „. dHj , „, , 
■(P)a + — ^(P)fe 



dx 



dy 



m 

dz 



{P)c 



dw 



(P)0 = deg{H,)Hi{P) = 0. 



Hence, 



^^\P),^{P),^{P) 



dHj 



■ u 



{P)a 



^^^(P)6+^(P)c: 



dy 



dz 



0. 



n 



dx dy dz J dx 

Thus u is orthogonal to all the planes IIj. 

Remark 5.2. Note that the previous lemma is also true for P singular and for each 
simple tangent not included at the plane at infinity. 

Applying the previous lemmas we get the next theorem. 

Theorem 5.3. Let Ex.E^ C C^ he such that 

1. El 

2. card(£'i°°) = card(^. 

Then, H(£^i,£^2) < oo. 

Proof. By (2), all points at infinity of both curves are simple and none tangent line is 
included at the plane infinity. Therefore, by Lemma 15. Ij all branches of the curves go 
to infinity following asymptotes. Since the direction vectors of the asymptotes depend 
only on the points at infinity (see Lemma [5. ip . by (1) the asymptotes of Ei and E2 are 
parallel. Now, the result follows reasoning as in the proof of Lemma 6.1. in [20]. n 



coo 
^2 



T) 



deg(^i) = deg(£2 



Remark 5.4. Note that, under the assumptions of Theorem 15.31 it also holds that 
H(£^i nlR^, £^2 nM^) < 00. We also observe that the previous theorem can be used under 
weaker assumptions for bounding H(<£^i fl M^,(£^2 H M^), since the result would follow 
analyzing the real asymptotes. 

Corollary 5.5. Let C,C be the input and output curves of our algorithm. Then 
H(C,C) < cx) anc^H(CnM^CnM3) < cx). 

Proof. It follows from the general assumptions, and from Theorems 13.11 and 15.31 n 
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6 Empirical Analysis of the Error 

We proved in the previous section that the Hausdorff distance H(C, C) is finite. In this 
section, because of the computational difficulties, instead of computing or bounding 
theoretically H(C,C) we approach the study of the Hausdorff distance empirically and 
apply it to the examples in Section HI Recall that 

H(C,C) := max{suppec{d(P,C)},supQ,c{d(g,C)}}, 

where C, C are the input and output curves of our algorithm. 

We explain next how to estimate suppg(^{d(P, C)} and supQgj{d((5,C)} indepen- 
dently. Let V(t) = {pi{t) ,p2{t) ,p3{t)) G R{tY be the parametrization of C output by 
our algorithm. 

Estimation of supQg^{d((5,C)} 

For every to ^ Q, such that V{to) is well deffned, and pr'(^o)p2'(^o)P3 '(^o) 7^ 0, we 
consider the normal plane to C at the point V{to) given by the parametrization: 

£i(to, h, fcs) = Vito) + kiVi{to) + k2V2{to), 

where fi(to) and f2(to) are unitary vectors in the direction of {—p3'{to),0,pi'{to)) and 
(0, — p3'(to),P2'(^o)) respectively. Moreover, we introduce the polynomials 

Viito, ki, /ca) = Fi(£i(to, ki, /ca)) E Q[ki, /^s], i = 1, . . . , s; 

we recall that {Fi, . . . , F,} are the polynomials generating C. In this situation, it holds 
that 

d(P(to), C) < min{||/cit;i(to) + A;2i'2(to)|| | A(to, k^, /ca) = 0, /ci, A;2 e C, z = 1, . . . , s}. 

Let pi(to) denote the r.h.s. of the previous inequality. 

We explain next how to choose an appropriate finite set T C Q to give 
maX(og7-{pi(to)} as an estimation of supQg^{d(Q,C)}. First, we obtain a finite set 
To C Q as follows. For each real pole of the parametrization V{t) we consider a finite 
sequence of isolating intervals {Jj}j=i,...,eo of length 1/10'^*+^'', and we take the middle 
point. Then To is the set containing all the middle points. Secondly we consider the 
set Ti = {(—2)* I 2 = 1, ... , 10*^1} for some ei G N and finally we take T = To U Ti. 

Estimation of suppgc{d(-P, C)} 

For every non singular point Pq = (a, b, c) G C, we consider the normal plane to C at 

C2{Po,x,y,z) = N,{Po){x - a) + N2{Po){y -b) + N,{Po){z - c) 



where (A''i(Po)) ^2(-Po)5 ^3(-Po)) is a unitary normal vector to C at Pq. Moreover let 
Q{Po,t) be the numerator of C2{Po,Pi(t):P2(t)^P3(t))- Then it holds that 

d(Po,C) < mm{\\Po~Vit)\\/g{Po,t) = and t G C}. 

Let P2(-Po) denote the r.h.s. of the previous inequality. 

We explain next how to choose an appropriate finite set £^ C C to give 
maxp(,g£-{p2(-Po)} as an estimation of suppg(.{d(P, C)}. For each i = 1, ... ,62, where 
62 G N, we compute the sets 00^, cof and u^ of intersection points of C with the planes 
X = (-2)\ y = {-2y and z = (-2)\ respectively. Then S = U-ii(w,^ U wf U cuf). It 
should be noticed that, in practice, the points in S cannot be taken as exact points on 
the curve but as e-points (see [I9]), thus contributing to an increment of the magnitude 
of the estimation. Observe that this is avoided when taking T. 

Distance Estimation for Examples in Section |4] 

Example 14. 2L The parametrization V{t) of C has 2 real poles. For eo = 20 we get 
maXio67-o{pi(^o)} = 0.2203928911 and for Ti = {(-2)^ | i = 1, . . . , 10^} we obtain 
maxjQg7-j{pi(to)} = 0.1036637452. Then our estimation of s\ipQ^(;{d{Q,C)} is equal to 
0.2203928911. 

For S = Uf£^(ci;f U wf U wf ) we obtain an estimation of suppg^l^l-P, C)} equal to 
maxp„gf {p2(i^o)} = 0.4705723389. Thus our estimation of H(C,C) is 0.4705723389. 

Example 14. 3L The parametrization V{t) of C has 2 real poles. For cq = 20 we get 
maxjoero{pi(^o)} = 0.1558549452 and for Ti = {(-2)' | i = 1, . . . , 10^} we obtain 
maxjQg7-i{pi(^o)} = 0.1603882181. Then our estimation of supQgj{d(Q,C)} is equal to 
0.1603882181. 

For S = Uf£^(a;f U wf U u^) we obtain an estimation of suppgc{d(P, C)} equal to 
maxp,g^{p2(Po)} = 0.1562381230. Thus our estimation of H(C,C) is 0.1603882181. 
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